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University  of  Wisconsin 
Madison,  Wisconsin  53706 

Near  the  end  of  the  nineteenth  century,  Boltzmann  (1894)  noted  in 
his  study  of  the  linear  diffusion  equation,  that  the  two  independent  vari- 
ables space  x and  time  t could  be  combined  into  a new  independent 
variable  | where  £ = £(x,  t)  . With  this  relation,  the  diffusion  equation 
which  is  a partial  differential  equation  (PDE)  could  be  transformed  into  an 
ordinary  differential  equation  (ODE).  Boltzmann  had  the  "Ansatz"  that 
the  new  variable  should  be  £,  = x/\/T  . 

At  approximately  the  same  time,  Sophus  Lie  (1881)  attempted  to 
construct  a general  integration  theory  for  differential  equations  utilizing 
the  theory  of  algebra.  Using  the  idea  of  continuous  groups  of  transfor- 
mations, he  was  able  to  reduce  the  order  of  an  ODE  and  in  some  cases 
obtain  a solution.  Although  not  developing  a general  integration  theory, 
he  also  examined  some  first  and  second  order  PDE. 

In  1952,  Morgan  (1952)  and  Michal  (1952)  presented  an  interesting 
simplification  of  Lie's  work  to  obtain  a complete  picture  of  the  group 
structure  of  partial  differential  equations.  Using  their  ideas,  as 
t On  leave  from  the  University  of  Iowa. 
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summarized  in  the  books  by  Ovsjannikov  (1962),  Ames  (1965,  1972)  and 
Bluman  and  Cole  (1971),  it  is  now  possible  to  remove  the  mystique  sur- 
rounding Boltzmann's  "Ansatz"  and,  in  fact,  have  a methodical  procedure 
to  construct  self  similar  variables  £ = £(x,  t)  and  to  examine  the  self- 
similar behavior  of  partial  differential  equations. 

An  alternative  approach  to  examine  the  self  similar  behavior  o 
PDE  is  to  use  the  ideas  of  dimensional  analysis  to  construct  the  seD 
similar  variables.  This  approach,  as  outlined  in  Sedov  (1959),  does  not, 
however,  yield  the  extensive  opportunities  that  the  group  theory  approach 
affords  and  will  not  be  discussed  further. 

To  date,  the  major  application  of  self  similar  analysis  has  been 
restricted  to  the  areas  of  fluid  mechanics  and  heat  transfer  where  exten- 
sive application  to  boundary  layer  phenomena  has  been  made.  Applica- 
tion to  physiological  and  electrical  circuit  problems  has  also 
recently  been  discussed  (Shen  (1976),  Lonngren  et  al.  (1975)).  Plasma 
physics  is  rich  in  phenomena  which  fall  into  the  self-similar  category  and 
should  receive  equivalent  attention.  To  date,  it  has  not.  It  is  to  this  end 
of  stimulating  our  colleagues  that  this  review  is  directed. 

In  Section  II,  the  detailed  procedure  to  obtain  the  self-similar  vari- 
ables and  the  resulting  ODE  will  be  presented.  We  shall  focus  on  the 
linear  diffusion  equation  as  being  the  vehicle  to  lead  to  an  understanding 
of  the  procedure.  Certain  constraints  imposed  by  either  boundary  or  initial 
conditions  or  by  conservation  laws  will  be  discussed. 


1 
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In  Section  III,  this  procedure  shall  be  applied  to  several  problems 
which  are  found  in  plasma  physics  and  that  we  have  treated,  mainly  with 
colleagues  at  Iowa.  They  include; 

(1)  One-dimensional  diffusion  where  the  diffusion  coefficient  is 
nonlinear  or  inhomogeneous  as  might  be  found  in  an  electron  cloud  ex- 
pansion, multipole  experiments,  or  in  studying  the  evolution  of  the  distri- 
bution function  of  particles  in  plasmas  supporting  various  aperiodic  in- 
stabilities; 

(2)  Diffusion  in  more  than  one  dimension  with  the  inclusion  of  a 
preferred  direction  of  drift; 

(3)  The  Korteweg-de  Vries  equation  which  can  be  used  to  model  ion 
acoustic  waves; 

(4)  Sets  of  ion  acoustic  wave  equations; 

(5)  The  Vlasov  model  of  a plasma. 

In  each  of  the  problems  that  is  discussed,  references  to  other 
theoretical  and  experimental  work  will  be  made. 

Self-similarity  is  not  the  panacea  to  solve  all  problems.  Some 
difficulties  that  we've  encountered  are  that  the  ODE  with  transformed 
boundary  conditions  may  not  be  amenable  to  solution,  neither  analytical 
nor  numerical.  Second,  if  solvable  mathematically,  the  solution  may  not 
describe  a physically  interesting  phenomena.  Third,  the  technique  is 
limited  to  problems  where  no  scale  length  nor  time  scale  such  as  fixed 
boundaries  exist  in  the  problem. 
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II  Self-Similar  Procedure 


In  order  to  find  the  similarity  variables,  we  make  use  of  a theory 
growing  from  the  Lie  theory  of  groups  where  it  has  been  shown  that  the 
similarity  variables  are  identical  to  the  invariants  of  a particular  one 
(or  more)  parameter  group  of  transformations.  We  shall  briefly  outline 
the  procedure,  details  and  references  can  be  found  in  the  texts  by  Ames 
(1965,  1972).  We  shall  examine  the  one-dimensional  linear  diffusion 
equation; 


P - p.  = 0 
xx  t 


(1) 


where  the  subscripts  denote  differentiation  with  respect  to  x and  t . 

We  shall  define  a one  parameter  ("a",  a is  positive  and  real)  group 

G as: 


G= , 


r Q - 

r P = a p 


x = aP  x 


t = a^t  . 


(2) 


This  is  called  the  "linear"  group.  Other  groups  exist  and  there  is  a 
"most  general"  group  called  the  "infinitesimal  group".  This  latter  group 
yields  all  possible  similarity  variables  but  one  is  soon  lost  in  a sea  of  al- 
gebraic formulas  that  it  is  difficult  to  discern  the  procedure  in  the  reams 
of  scratch  paper.  One  can  refer  to  Ames  (1965,  1972),  Bluman  and  Cole 
(1974)  or  Shen  (1976)  for  examples.  It  is  not,  however,  entirely  hopeless, 
as  computer  programs  can  now  be  written  to  handle  large  arrays  of 


-4- 


simultaneous  symbolic  algebraic  equations. 


In  (2),  a,  [3,  and  y are  constants  which  are  determined  such  that 
(1)  is  "(absolutely)  constant  conformally  invariant"  under  the  group  G 
(Ames  196s,  1972).  A function  F(y)  is  said  to  be  "constant  conformally 
invariant"  (CCI)  under  G if  F(y)  = f(a)  F(y)  where  f(a)  is  some 
function  of  the  parameter  a . If  f( a)  h 1,  the  constant  conformal  in- 
variance is  called  "absolute".  (ACCI) 

Substituting  (2)  in  (1),  we  write 


a- 23  - «-Y  - „ 

a rp-  - - a p-=0. 
xx  t 


For  (3)  to  be  ACCI  under  the  transformation  group  G one  requires 

a - 20  = a - y (4) 


or  y = 2(3  . We  shall  defer  until  later  the  further  specification  of  these 
constants. 

Instead,  we  now  seek  to  determine  the  "invariants"  of  the  trans- 
formation group  G . This  is  achieved  by  employing  a theorem  from  group 
theory  (Ames  (1965),  (1972)).  The  invariants  are  obtained  from  Ql  = 0 
where  I is  an  invariant  and  0 is  the  operator 


0- 


a 

9x 

8 

at 

a=l  9p 

+ 8a 

8x 

a=l 

+ 8T 

d d d 

-“p  57  - Tx  - vt  a?  ' 


The  solutions  of  QI  ^ 0 are  obtained  by  solving  the  Lagrange  subsidiary 
equations 

dp  _ _dx_  _ _dt_ 

-ap  -px  -yt 

According  to  the  theorem  developed  by  Morgan  (1952),  these  "in- 
variants" are  the  self  similar  variables.  Solutions  of  (6)  are 


<M£)  = 


P(x,  t) 
t °/y 


and 


(7) 


From  (4),  we  found  that  p/y  = 1/2  . Therefore  the  Boltzmann  transformation 
is  recovered.  Note  that  we  could  have  combined  (6)  in  a different  order 
and  obtained 


^ and  = 


(8) 


Having  found  the  self-similar  variables,  let  us  transform  the  PDE 
(1)  using  the  self  similar  variables  (7)  into  an  ODE.  The  result  is 


or 


4>  + f 4>  - - 4>  = ° 

££  2 £ V 


+2(£/Z)  4»e/2  -4  - <t>  = 0 

2 2 


(9) 


(10) 


Writing  the  ODE  in  the  form  (10)  allows  us  to  recognize  that  its  solution 
can  be  written  in  terms  of  complementary  error  functions  (Gautschi  (1964)) 


2- 


Qt 

2- 


4>  = Ai  ^ erfc(  |-  ) + Bi  ^ erfc(-  |/2) 


(11) 
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where  i 


is  an  ordering  parameter  and 


.-1  , t /,  2 -4  /4 

1 erfc  4/2  = c 


si  TT 


i°  erfc  4/2  = erfc 


(12) 


n n 1 

in  erfc  4/2  = f in_1  erfc  tdt  n = 0,  1,  2, 

4/2 


At  this  stage,  the  parameter  a/y  is  still  arbitrary.  We  shall  specify 
it  to  satisfy  boundary  conditions  or  a conservation  law.  We  note  that  two 
boundary  conditions  on  p(x,  t)  have  necessarily  "consolidated"  into  one 
for  6,  namely 


p (x  = »,  t)  = 0 

P(X,  t = 0)  = 0 


4>(4=x)  = o 


(13) 


The  third  boundary  condition  could  have  one  of  two  forms  which  would  yield 
self  similar  solutions.  They  are 

p(x  = 0,  t)  = constant  (i4) 

00 

J pdx  = constant  . (15) 

0 


Since  p(x  = 0,  t)  transforms  via  (7)  to  b(e  = 0)>  we  note  that  (14)  requires 
that  a/y  - 0 . The  self  similar  solution  (11)  for  this  boundary  condition 
is 


-7- 


” 


p = <?  A erfc(£/2)  A erfc(  — ^ ) (16) 

zn/t 

where  the  constant  B = 0 in  order  to  satisfy  (13). 

Although  not  proven  to  our  knowledge  but  generally  accepted  as  a 
"Similarity  Postulate"  (Moran  and  Gaggioli  (1969)),  the  conservation  law 
given  in  (15)  should  be  invariant  under  the  group  transformation  in  order 
to  have  similarity  solutions.  Applying  (2)  to  (15),  we  write 

or  oc 

J pdx  = aa  ^ J p dx  . (17) 

0 0 

For  this  to  be  ACCI,  we  have 


a £ _ i 
Y ~ Y 2 

The  self  similar  solution  (11)  that  satisfies  this  conservation  law  and  (13) 
is 


<(>  2 A'  -x  /4t 

p(X,  t)  = — = _ ■_  « 

\I  t \l  IT  sl  t 


(18) 


where  the  constant  A'  can  be  determined  by  evaluating  (15)  using  (18). 

The  procedure,  as  described  in  this  section  will  be  applied  to 
several  examples  found  in  plasma  physics  in  the  next  section. 


-8- 


Ill  Applications  to  Plasma  Physics 

Plasmas  are  a particularly  interesting  and  fruitful  area  in  which  to 
:pply  the  self  similar  technique.  This  is  particularly  true  when  it  can  be 
argued  that  there  is  no  scale  length  that  is  important.  However,  there 
usually  is  a relevant  scale  length  such  as  the  Debye  length  in  most  plasma 
problems,  at  least  in  some  initial  stage.  The  se!f  similar  analysis  may, 
in  those  cases,  yield  valuable  information  for  some  asymptotic  state 
(Barrenblatt  and  Zel'dovich  (1972)). 

In  the  following,  we  shall  select  several  examples  from  plasmas 
which  fall  into  the  self-similar  class.  To  analyze  each  problem,  we  follow 
the  procedure  given  in  Section  II  and  only  describe  the  physical  phenomena, 
list  the  PDE,  the  self  similar  variables,  the  ODE,  and  if  possible,  the 
solution  to  the  ODE  without  repeating  the  procedure  for  each  problem. 
References  for  each  problem  are  also  presented. 

A.  Linear  1-D  Diffusion  Equation,  Inhomogeneous  Diffusion  Coefficient 

In  studies  of  the  evolution  of  the  distribution  function  of  particles 
in  plasmas  supporting  various  aperiodic  instabilities  or  in  calculations 
involving  the  forward  scattering  of  photons  by  plasmons,  it  has  been  found 
that  the  problem  could  be  modeled  with  a diffusion  equation  (Peyraud  and 
Coste  (1976)). 

Pt  = ^ Px  lx  • (19) 

Using  the  linear  group  G (2),  we  find  the  self  similar  variables  to  be 


¥ 
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Equation  (19)  transforms  to 


and  d = 


- $ + — =7  4 »•  - (£  m 4>  ) • 

V m-2  4 4 s 


The  requirement  of  "consolidation"  (13)  specifies  that  = — - — > 0 or 

•y  2 -m 

m < 2 . 


The  constant  a/y  can  be  specified  to  satisfy  two  boundary  con- 
ditions, namely  (14)  and  (15).  The  solutions  for  the  two  cases  p(0,t) 

30 

constant  and  J pdx  = constant  are  (Lonngren  ( 1 9 7 6 ) > ; 

0 


p(x,  t)  = K 


[-(S)J 


de  - 1 


k r x2-m 

^ exp  T 

t L.  t(2-m) 


respectively  and  where  K is  a constant. 

As  the  diffusion  equation  can  be  also  used  to  model  a distributed 
RC  transmission  line  (eg.  Lonngren,  et  al.  (1975)),  there  can  be  a practical 
application  of  this  result.  In  particular,  let  us  assume  that  all  the  capacitors 
have  the  same  value  and  the  resistors  are  distributed  inhomogeneously 
such  that  R ~ RQ  x m . If  such  a line  were  constructed  and  it  were  suf- 
ficiently short  such  that  exponential  term  in  (2  3)  remained  approximately 


constant,  then  the  voltage  response  of  the  line  to  an  impulse  source  would 
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be  proportional  to  1/t^  " . Such  a line  has  been  used  in  biomedical 

applications  for  the  value  m = 1 (Gagne  and  Poussart  (1976)). 

B.  Nonlinear  1 - D Diffusion  liquation 

Several  phenomena  in  plasma  physics  can  be  modeled  with  the 
nonlinear  diffusion  equation. 

[pn  Px  lx  ' Pt  = 0 • (Z4) 

For  example:  (1)  Recent  experiments  in  multipoles  have  confirmed  that  in 

certain  regions,  the  diffusion  coefficient  for  particles  across  the  magnetic 

-1/2 

field  depends  on  (the  density  of  the  particles)  . (Drake  (1973), 

Berryman  (1976)); (2)  In  studying  the  skin  current  penetration  into  turbulent 
plasmas  where  the  conductivity  depends  on  (the  local  electric  field)1"1 
(Hirose  et  al.  (1970)),  it  has  been  shown  that  the  governing  equation  is 
of  the  form  of  (24)  (Hirose  and  Alexeff  (1973)).  An  examination  of  this 
equation  allowed  us  to  comment  on  penetration  times  under  various  stages 
of  ion  acoustic  turbulence  (Lonngren,  et  al.  (1974)).  (3)  As  the  conductivity 
of  the  turbulent  plasma  saturated  with  increasing  field  (Hirose  et  al.  (1970)), 
we  could  model  the  conductivity  as  & ~ tr  exp(-E)  and  have  been  able 
to  obtain  (24)  with  n = -1  . (Anmadi,  et  al.  (1976)).  A similar  equation 
describes  the  expansion  of  a Maxwellianized  electron  cloud  into  a 
vacuum  (Lonngren  and  Hirose  (1976)). 

For  these  problems,  the  self  similar  variables  are  of  the  form 
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t = 7^  and  * = 


n (2  -1) 

n y 


where  the  parameter  p/y  is  chosen  to  satisfy  the  boundary  conditions  or 


conservation  laws  and  4>  satisfies 


(4>V),  + £ £«>  - - (2&  - l)<t*  = 0 

4 4 Y 4 n y 


Again  the  boundary  conditions  (14)  and  (15)  can  be  applied.  Using 
(15),  we  find  that  p/y  = - ^ (2  p/y  - 1)  and  (26)  can  be  directly  integrated 
(Ames  (1965),  Gilding  and  Peletier  (1976))  to  yield 


H>] 


2>>  - 
1 n 


i < 4, 


where  £ 


|~2(n+2)  1//2 


n > 0 . 


This  could  be  considered  a "sharpfront"  solution  in  that 


*n+1(4  = 40)  = o 


For  -2  < n < 0,  say  n = -1;  the  solution  is 
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<?  = 


1 

k +tz/i 


or  p = 


1 


r 2 

1 ^tK 
Lkt 


(29) 


1 Q 

which  is  valid  for  - ^ - 0 as  x - <«  (Lonnqren  and  Hirose  (1976)). 

p dx 

To  apply  the  second  boundary  condition  (14),  it  is  convenient  to 
let  d - i ' ''  and  note  that  the  constant  (3/y  - 1/2  . Equation  (26)  be- 
comes 


-Hi.. 


(£  ' £ ' 0^  + J 6 


= 0 


(30) 


The  solution  for  n = -1  with  the  boundary  conditions  d(|  = 0)  = 1 and 
O(|=ao)  = 0 has  been  treated  by  Fujita  and  is  described  in  the  classic 
book  by  Crank  (1964).  Approximate  solutions  were  recently  discussed  by 
Ahmadi,  et  al.  (197  6). 

Solutions  of  other  one-dimensional  diffusion  equations  have  been 
given  by  e.g.  Phillip  (I960),  Singh  (1967),  Swan  (1976)  and  Tuck  (1976) 
and  should  be  referred  to. 


C.  Diffusion  In  More  Than  One  Dimension 

In  several  cases  in  plasma  physics,  it  is  important  to  examine 
diffusion  of  particles  in  more  than  one  dimension  and  in  cases  where  the 
parallel  and  perpendicular  diffusion  coefficients  are  not  equal,  parallel 
could  be  interpreted  as  being  along  an  external  magnetic  field.  In  such 


cases,  the  linear  diffusion  equation  becomes 


~ 2 


Pt  + WPX=  Dn  Pxx  + Dj.  V 


(31) 
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r 


where  w is  a drift  velocity.  As  and  are  assumed  at  this  stage 

to  be  linear,  we  can  renormalize  the  transverse  coordinates  such  that 


Di^ip=  Dn  vlp  • 


(32) 


In  addition,  letting 


( 31)  becomes 


•& 


w w 

P=  N exp|  — x - ^5—  t 
U11  *U11 


T = Dllt 


■1 


(33) 


N = N + V N 
t xx  ± 


(34) 


Using  the  procedure  outlined  in  Section  II  with  the  addition  of  two 
terms  to  the  group  G (2),  it  follows  that  the  group  invariants  are 


. x , _ y z . _ N 

i - > 4 - > ^ - / 

\]  T \/  T \/  T T ^ 


(35) 


where  y and  z are  in  the  transverse  direction.  Substituting  (3  5)  in 
(34),  we  obtain 


^ «>  - “ fe4>e  + = + V±<)> 


(36) 


where  V 2 indicates  V2  in  the  transformed  transverse  variables  £,  and 
± X 


We  shall  examine  the  conservation  law  only  (15).  For  this  problem, 
it  generalizes  to 
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J N dxdydz  = constant  (37 

which  implies  that  a/y  - - 3/2  . Using  this  in  (36)  suggests  that  we 


j[U4>)e  + (M>);  + ] = 4>e.  + V^4» 


which  is  amenable  to  a treatment  using  separation  of  variables  and  inte- 
gration. Note  that  if  the  problem  were  posed  in  cylindrical  coordinates, 
the  proper  conservation  law  would  be 


J Ndrdx  = constant  (39) 

as  there  is  a fixed  scale  length  of  2ir  in  the  third  coordinate.  This  im- 
plies that  a/y  = -1  . 

The  procedure  that  follows  is  straight  forward  and  will  not  be  re- 
produced here.  In  cylindrical  coordinates,  we  find 


const 

(Dllt) 


t 2 \ r-  2 2 - 

T0  \ 8D  t ) eXP[_'  8D  t ‘ 4Djjt  + 2Dn(X  2 J 


This  agrees  with  a result  of  Eastlund  (1966),  who  examined  diffusion  in 
a O-machine. 

This  diffusion  equation  (31)  has  found  importance  in  several  non- 
plasma areas  also.  We  cite  as  examples:  Dispersion  in  uniform  porous 
media  flow  of  ground  water  (Shen  (1976))  and  in  the  spread  of  cancer  in 
the  uterus  (Swan  (197  5)).  Extensions  to  nonlinear  and  inhomogeneous 
problems  would  be  worthwhile  as  we've  observed  in  the  one  dimensional 


case. 
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D.  Korteweg  de  Vries  Equation 


Considerable  attention  has  been  given  to  understanding  the  Korteweg- 
de  Vries  equation. 

pt  f ppx  f pxxx  = ° • (41) 

For  example,  Washimi  and  Taniuti  (1966)  showed  that  the  low  frequency 
ion  acoustic  wave  could  be  modeled  with  this  equation  as  a first  approx- 
imation for  including  nonlinear  effects. 

Using  the  linear  group  G,  given  in  (2),  one  computes  the  self  similar 
variables  to  be 

s=773  ond  ♦'“dyr  (4Z) 

where  4>  satisfies 

-|  4>  - y £ + <M>£  + 6 2 = 0 . (43) 

A discussion  of  the  solution  of  (4  3)  is  given  in  Berezin  and  Karpman  (1964). 

Two  extensions  to  this  solution  have  been  given.  For  cases  where 
the  nonlinear  term  pp^  in  (41)  can  be  neglected,  the  self  similar  solutions 
for  (41)  can  be  written  in  terms  of  Airy  functions  or  integrals  of  Airy 
functions  which  have  the  argument  of  the  similarity  variable  (42).  These 
two  solutions  satisfy  the  conservation  law  (15)  and  the  boundary  condition 
(14)  respectively  (Ikezi,  et  al.  (1974)). 

The  second  extension  by  Shen  and  Ames  (1974)  uses  the  infinitesimal 
group  rather  than  the  linear  group  that  we've  used  and  has  ascertained  all 
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possible  similarity  variables  for  the  KdV  equation.  They  find 
<t>  = (at  + 6)*”  (u  - — ap)  ^ a ^ 

(44) 

^ = 7 a - 3u(at+6)  + b(  ap-a&)]  /(at+6)  * 3 

where  a,  a,  p,  and  6 are  constants  and  4>  satisfies  (43) 

E.  Ion  Acoustic  Wave  Equations,  Fluid  Model 

In  a series  of  two  papers  (Hsuan,  et  al.  (1974)  and  Shen  and  Lonngren 
(1976a)),  the  similarity  properties  of  four  sets  of  fluid  equations  that  have 
been  used  to  describe  the  propagation  of  ion  acoustic  waves  in  a plasma 
are  discussed.  The  sets  of  equations  are;  (a)  multiple  species  fluid 
equations  truncated  at  the  third  moment  plus  Poisson's  equation; 

(b)  massless  isothermal  linear  electron  fluid  and  cold  nonlinear  ion  fluid 
plus  Poisson's  equation;  (c)  massless  isothermal  nonlinear  electron  fluid 
and  cold  nonlinear  ion  fluid  plus  Poisson's  equation;  and  (d)  massless 
isothermal  electron  fluid  and  cold  nonlinear  ion  fluid  with  a quasineutrality 
assumption. 

Set  (c)  is  the  most  general  and  the  infinitesimal  group  was  applied 
in  addition  to  the  linear  group  G (2).  (Shen  and  Lonngren  (1976a)  Shen 
(1976)).  As  a result  several  possible  similarity  variables  were  found  which 
included  those  quoted  by  Zhmudsky  (197  5)  and  the  traveling  wave  variables. 
Numerical  difficulties  prevented  a complete  integration  of  the  resulting  of 
ODE  (Shen  (1976)).  Also  a scale  length  (the  Debye  length)  existed  in  the 
problem. 


Set  (d)  has  been  frequently  studied  in  plasma  problems  (eg.  Alexeff 


et  al.  (1971),  Allen  and  Andrews  (1970)).  The  equations  are  identical  to 
those  in  ordinary  fluids. 

Finally,  set  (c)  was  examined  in  a nonneutral  plasma  approximation 
where  the  effect  of  the  electrons  was  neglected  (Shen  and  Lonngren  (1976b)). 
Such  a model  extends  a paper  by  Gintsburg  (1974)  and  is  germane  to  re- 
lativistic electron  beam  devices. 

Below,  we  tabulate  the  sets  of  PDE,  the  similarity  variables  and  the 
insulting  ODE.  The  results  for  sets  (a)  and  (d)  were  obtained  using  the 
"linear"  group,  (b)  was  obtained  using  the  "spiral"  group  and  (c)  was 
obtained  using  the  "infinitesimal"  group.  Details  appear  in  (Hsuan  et  al. 
(1974),  Shen  (1976)  and  Shen  and  Lonngren  (1976a)). 


a)  The  PDE  are: 


3n 


1 


at 

9v 


<57<Vii  = o , 


_i  t + j_  !!l  . q_±  E 

at  3x  n^  3x  m^ 

aPj  avl  aPj 

“aT  + Yipi  IbT  + vi  “ax'  = 0 * 


(45) 


£E_  _ v qi 

8x  - I .o"1  • 


The  similarity  variables  are; 


L 
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The  similarity  variables  are: 

£ = t exp(- ax),  N.  = n.  exp(2ax)  , 

N = n exp(2ax),  U = v exp(ax)  , 
e e i 


(49) 
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The  ODE  are: 


1 

dN. 

i 

dU 

e - eu) 

a 

d£ 

- £N.  TT 
l d£ 

1 

dU 

dN 

e 

(-  - £U) 

Or 

d£ 

‘ ^ d£ 

? 

, d N 

dN 

.2  e 

5 

+ 5£ 

H * " 

2 i 


(SO) 


c)  The  PDE  are: 


9n  9(n  ,v.) 

-1  + — 

9t  9x 


= 0 


9v . 


9v 


9t  + Vi  9x 


92<$> 

2 ~ ne  ' ni 


9x 


9<t> 

9x 


(51) 


n = e 
e 


The  similarity  variables  are; 


£ = 


X+t 

V 


t + * 


v 


ln(t  + -£■)  + 


2 

v 


t +J 


Ni<£)  = (t  + -£)  n.(x,t)  , 

Vi(£)  = v1(x,t)  - ~ In ( 1 +f)  , 


(52) 


Ne(£)  = (t  + ng(x,t)  . 
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The  ODE  are; 


-2N.  - 


(V.  -£) 


edN. 

d(N.  V.) 

l l 

6 

dN. 

l 

dE 

d£ 

V 

d£ 

dV. 

l 

, dN 

1 e 

dV. 

d£ 

f N d£ 

e 

V 

( de 

l w 1 ^ 

— (—  - ) = N - N 

d£  (Ne  d £ ’ e i 


d)  The  PDE  are-. 


— + = o , 

at  9x 


3v  dv  A 3n 

at  dx  nQ  3x 


The  similarity  variables  are: 


The  ODE  are; 


u = 7rr>  N=7iirr' 


-A4  + 2CA-DN  + = 0 , 


„ du  TT  dU  A dN 
(A-l)U  - A4  de  + u de  + ^ 


F.  Vlasov  Model  for  Ion  Acoustic  Waves 

The  distribution  function  for  ions  in  the  absense  of  a magnetic  field 
can  be  computed  from  the  Vlasov  equation 


3f  3f  e_  £f  _ q 
Jt  + V 3y  ” M 9y  3v 
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With  a quasi-neutrality  assumption,  n ~ n,  and  with  a Boltzmann  ap- 

e 1 

proximation  for  electrons  n = n„  expfe  o /T  1,  we  write  (57)  as 

e 0 L e 


9f  + u iL  IL  d0n  I fdU)  _ 
at  ax  ‘ au  dx 


(58) 


Standard  dimensionless  parameters  are  used  in  (58),  namely 

n ,,  y oo 

e TT  V p 

n = — , t = w . t U = , x = — 

0 V T_  s]  T 


e/M 


e/M, 


Applying  the  linear  group  G (2)  to  (58),  we  find  that  the  self 
similar  variables  are 

f 


| = - and  4>  = i~ 

t t«/Y 


(59) 


where  4>  satisfies 


a d^_  d^_  d4>_  d(ln  J <t>dU) 

Y ^ d£  d£  ' dU  d£  ’ 


(60) 


Again,  we  have  two  possible  boundary  conditions  that  can  be  treated  ; 
(14)  and  (15).  In  a pioneering  series  of  papers,  Gurevich  and  Pitaevsky 
examined  the  solutions  governed  by  boundary  condition  (14),  i.  e. 
f(x  = 0,  t)  = constant.  Their  work  has  been  summarized  in  a recent  review 
paper  (Gurevich  and  Pitaevsky  (197  5)).  Experiments  performed  by  Korn, 
Marshall  and  Schlesinger  (1970)  and  others  seem  to  confirm  the  predictions. 

The  second  boundary  condition  (15),  J f(x,  t)dx  = constant  has  re- 
ceived considerable  attention  from  Jensen  and  his  colleagues  in  their 
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investigation  of  the  Green's  function  for  a linearized  version  of  (58) 


(Anderson  et  al.  (1971)  Jensen  et  al.  (1974)  Christoffersen  et  al.  (1974)). 
They  found  that  the  leading  term  in  the  series  for  the  Green's  function  is 
the  self-similar  term  --  h(^-)  which  is  noted  by  setting  a/y  = -1  in  (59). 
Equation  (60)  in  this  case  is  written 


IV  Conclusion 


In  this  tutorial  and  review  paper,  we  have  summarized  the  technique 
of  self-similar  solution  of  partial  differential  equations  and  presented 
several  examples  where  it  has  been  applied  to  problems  found  in  plasma 
physics. 

The  author  wishes  to  acknowledge  H.  Hsuan,  H.  Shen,  A.  Hirose, 
and  W.  Ames  who  actively  contributed  to  this  work. 
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